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We  numerically  demonstrate  optical  negative  refraction  in  ferrofluids  containing  isotropic  Fe304 
nanoparticles,  each  having  an  isotropic  Ag  shell,  in  the  presence  of  an  external  dc  magnetic  field  H. 

The  all-angle  broadband  optical  negative  refraction  with  magnetocontrollability  arises  from  //-induced 
chains  or  columns.  They  result  in  hyperbolic  equifrequency  contour  for  transverse  magnetic  waves 
propagating  in  the  system.  The  finite  element  simulations  verify  the  analyses  using  the  effective  medium 
approximation.  Experimental  demonstration  and  potential  applications  are  suggested  and  discussed. 
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In  1968,  Veselago  [1]  theoretically  investigated  the  elec¬ 
trodynamic  consequences  of  a  medium  with  simultane¬ 
ously  negative  permittivity  and  permeability.  He  predicted 
that  such  a  medium  possesses  a  negative  (phase)  index, 
which  can  result  in  a  reversed  Snell’s  law,  i.e.,  negative 
refraction.  People  have  realized  negative  refraction  of 
optical  waves  or  microwaves  in  different  systems,  includ¬ 
ing  metamaterials  [2,3],  photonic  crystals  [4],  plasmonic 
waveguides  [5],  chiral  media  [6],  and  superconductor- 
ferromagnet  superlattices  [7].  The  significant  application 
of  negative  refraction  is  the  concept  of  a  perfect  lens, 
which  can  lead  to  subwavelength  imaging  beyond  the 
diffraction  limit  [8]. 

However,  almost  all  the  existing  methods  for  achieving 
negative  refraction  were  proposed  or  established  in  the 
realm  of  solid  materials,  in  contrast  to  soft  materials  with 
the  specific  characteristic  “softness”.  Literally  softness 
might  offer  an  extra  freedom  of  tailoring  physical  proper¬ 
ties;  hence,  it  encourages  us  to  investigate  optical  refrac¬ 
tion  in  certain  soft  materials.  As  a  result,  we  reveal,  for  the 
first  time,  a  new  class  of  all-angle  broadband  optical  nega¬ 
tive  refraction  in  ferrofluids  with  magnetocontrollability. 
Its  underlying  mechanism  arises  from  assembly  metallic 
chain  or  column  structures  induced  by  an  external  dc 
magnetic  field  H.  This  work  paves  a  new  way  for  designing 
tunable,  active  metamaterials. 

In  general,  ferrofluids  are  colloidal  suspensions  com¬ 
posed  of  ferromagnetic  or  ferrimagnetic  nanoparticles  of 
about  10  nm  diameter  dispersed  in  a  carrier  fluid,  usually 
water  or  kerosene  [9].  Two  typical  kinds  of  materials  for 
fabricating  such  nanoparticles  are  Co  and  Fe304  [9].  In  this 
work  we  shall  consider  an  aqueous  ferrofluid  system  which 
contains  spherical  isotropic  Fe304  nanoparticles,  each 
coated  with  a  spherical  isotropic  Ag  shell,  in  the  presence 
of  H  in  z  axis.  The  reason  why  we  choose  Fe304  rather 


than  Co  is  due  to  its  lower  optical  absorption.  In  principle, 
the  Ag  used  for  making  shells  may  be  replaced  by  other 
metallic  materials  having  a  permittivity  with  a  big  negative 
real  part  and  small  positive  imaginary  part.  If  water  is 
replaced  with  kerosene,  the  qualitative  results  as  to  be 
revealed  would  remain  unchanged.  Both  experiments 
[10-12]  and  computer  simulations  [13,14]  have  shown 
that  nanoparticle  chains  can  be  formed  in  ferrofluids  as 
the  magnetic  dipolar  interaction  overwhelms  the  thermal 
energy.  For  ferrofluids,  the  field-induced  columnar  phase 
with  equal  spacing  was  experimentally  reported  [11,12]  in 
confined  ferrofluid  films  subjected  to  an  in-plane  //,  and 
the  dipolar  interaction  between  the  chains  (which  are  lo¬ 
cally  displaced  in  a  hexagonal  fashion)  gives  rise  to  their 
arrangement  in  columns.  In  other  words,  nanoparticles  in 
ferrofluids  can  exist  in  the  form  of  both  (either)  chains  and 
(or)  columns.  In  our  work,  the  word  “chain”  denotes  a 
single-particle- width  z-directed  array  of  touching  Fe304 
nanoparticles  with  an  Ag  shell,  and  “ column ”  represents 
the  aggregation  of  such  chains,  namely,  the  chains  aggre¬ 
gate  together  to  form  z-directed  thick  columns. 
Incidentally,  throughout  this  work  both  permittivities  and 
permeabilities  denote  relative  ones. 

An  effective  medium  approximation  and  a  2D  finite 
element  method . — Let  us  start  by  considering  the  ferrofluid 
system  containing  Fe304  nanoparticles  coated  by  an  Ag 
shell,  in  the  presence  of  z-directed  H.  The  chains,  each 
having  the  same  number  n  of  nanoparticles  [15],  are  in¬ 
duced  to  appear  by  H.  They  are  all  directed  along  z  axis, 
and  are  assumed  to  be  randomly  suspended  in  the  system. 
In  the  quasistatic  approximation,  we  can  utilize  the  aniso¬ 
tropic  form  of  the  effective  medium  approximation  (EMA) 
[16]  to  calculate  the  effective  permittivity  tensor  with  three 
nonzero  diagonal  components  sxx ,  sxx),  and  szz 
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h  +  (v~  l)sx 


+  (1  ~p) 


^2 


e2  +  -  l)sx 


=  0, 


(1) 


where  p  is  the  volume  fraction  of  the  coated  nanoparticles, 
and  8  2  the  permittivity  of  water.  Except  for  the  EM  A,  some 
other  theories  which  appear  to  be  more  complex  can  also 
be  adopted,  e.g.,  the  spectral  representation  approach 
[16,17].  For  Eq.  (1)  we  have  modeled  an  individual  chain 
as  a  spheroid  with  gx  (shape  factor  perpendicular  to  H)  and 
gz  (shape  factor  along  H)  [18].  They  satisfy  a  geometrical 
sum  rule  gz  +  2 gx  =  1  [19].  It  is  important  to  note  that  the 
result  calculated  by  the  isotropic  counterpart  of  the  EMA, 
namely,  Eq.  (1)  with  gxz  =  1/3,  was  demonstrated  to 
reasonably  lie  in  the  Hashin-Shtrikman  bounds  which  pro¬ 
vide  the  tightest  constraints  for  the  effective  permittivity  of 
a  composite  comprising  an  isotropic  mixture  of  two  iso¬ 
tropic  dielectric  materials  [20].  Thus  in  view  of  the  physi¬ 
cal  reasonability  and  formal  compactness,  we  prefer  to  use 
the  EMA  [Eq.  (1)]  for  treating  the  ferrofluid  system.  In 
Eq.  (1)  si  denotes  the  equivalent  permittivity  of  the  coated 
nanoparticle.  It  can  be  given  by  solving  Laplace’s  equation 
of  electrostatics  together  with  appropriate  boundary  con¬ 
ditions  [21],  which  is  valid  in  our  quasistatic  system  where 
the  incident  wavelength  A  is  much  larger  than  the  size  of 
the  coated  nanoparticles.  After  some  straightforward  deri¬ 
vations,  we  obtain 


8/(1  +  2 p)  +  2sa(l  -  p) 
£/(l  -  p)  +  £a( 2  +  p) 


(2) 


with  p  =  r3/(r  +  cl)3,  where  r  is  the  radius  of  the  core 
(Fe304)  with  permittivity  Sf,  and  d  is  the  thickness  of  the 
shell  (Ag)  with  permittivity  sa.  Equation  (2)  is  an  exact 
solution  for  an  isolated  Fe304  nanoparticle  with  an  Ag 
shell,  but  it  serves  as  the  result  of  the  first-order  approxi¬ 
mation  for  such  a  nanoparticle  in  our  ferrofluid  system  by 
neglecting  local  field  effects  caused  by  all  the  other  nano¬ 
particles.  In  Eq.  (1),  gz  is  given  by  a  simple  approximation 
gz  =  1/(1  —  n 2)  +  n\n(n  +  Vn2  —  1  )/(n2  —  l)3/2  [22], 

where  n  >  2.  In  reality,  n  is  proportional  to  the  ratio  of 
magnetic  energy  to  thermal  energy  (a  lowest-order  as¬ 
sumption).  Thus,  without  loss  of  generality,  n  will  be 
used  to  equivalently  represent  the  strength  of  H.  This 
correlation  corresponds  to  the  fact  that  higher  H  yields 
larger  n  [13].  Apparently  H  =  0  causes  n  =  1  or  gz  = 
1/3.  In  this  case,  the  nanoparticles  are  randomly  suspended 
in  the  ferrofluid  system,  resulting  in  an  isotropic  effective 
permittivity.  It  is  worth  noting  that  the  EMA  [Eq.  (1)]  has 
included  energy  dissipation  related  to  absorption.  On  the 
other  hand,  the  effect  of  scattering  is  assumed  to  be  ne¬ 
glected  due  to  the  small  scattering  cross  section  of  a  chain 
with  respect  to  an  incident  wavelength. 

We  use  the  software  COMSOL  Multiphysics  3.5  to  per¬ 
form  finite  element  simulations  of  optical  refraction  in 
ferrofluids,  see  Fig.  1.  A  transverse  magnetic  (TM) 
Gaussian  beam  with  a  width  of  1.6  fi m  in  vz  plane  with 
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FIG.  1  (color  online).  2D  finite  element  simulations  of  the 
distribution  of  the  absolute  value  (denoted  by  the  color  or 
brightness)  of  electric  field  components  of  an  incident  TM 
wave  with  A  =  758  nm,  in  a  ferrofluid  system  of  thickness 
1  yum  containing  chains  of  Fe304  nanoparticles  coated  by  an 
Ag  shell.  The  ferrofluid  system  is  replaced  by  an  effective 
medium  described  by  the  parameters  calculated  according  to 
Eq.  (1).  The  direction  and  size  of  the  blue  or  dark  arrows 
indicate,  respectively,  the  direction  and  magnitude  of  the  local 
power  flow  or  Poynting  vector. 


a  y-directed  magnetic  field  component  is  incident  to  the 
ferrofluid  system  at  an  angle  <f>  =  30°  with  respect  to  z 
axis.  Here  </>  just  denotes  the  angle  of  incidence.  Using 
Eq.  (1),  we  obtain  sxx  =  4.43  +  0.15/  and  szz  =  —3.29  + 
0.19/  according  to  the  parameters:  =  —21.4  +  0.9/, 

s2  =  1.77,  p  =  0.22,  n  =  50,  and  the  length  of  each  chain 
is  1  pum.  Here  has  been  calculated  according  to  Eq.  (2) 
with  the  parameters:  r  =  5  nm,  d  =  5  nm,  Sf  = 
4.8  +  2.2/  (as  experimentally  measured  at  A  =  758  nm 
for  magnetite,  a  ferrimagnet  and  crystallites  with  the  in¬ 
verse  cubic  spinel  structure  [23]),  and  sa  =  —20.61  + 
1.27/  (which  was  measured  for  a  Ag  thin  film  at  A  = 
758  nm  [24]).  Note  we  already  take  into  account  the  con¬ 
ductivity  property  of  the  materials  by  using  the  complex 
permittivity  values  since  the  conductivity  and  permittivity 
are  related  [25].  Figure  1  plots  the  distribution  of  the 
absolute  value  of  electric  field  components,  and  it  clearly 
demonstrates  the  appearance  of  negative  refraction. 

To  understand  the  finite  element  simulation  results  de¬ 
picted  in  Fig.  1,  we  resort  to  a  theoretical  approach  for  the 
system  with  Re(e  «)<o  and  Re(8xx  xy)  >  0.  Here  Re(-  •  •) 
means  the  real  part  of  •  •  • .  On  the  basis  of  Maxwell’s 
equations,  for  the  TM  waves  with  the  magnetic  field  com¬ 
ponent  polarized  in  the  y  axis  and  the  electric  field  com¬ 
ponent  located  in  the  vz  plane,  the  dispersion  relation  for 
the  wave  propagating  in  a  general  anisotropic  medium  is 

(k2Jezz)  +  (k2z/ sxx)  =  P'yyko,  (3) 

where  kx  (or  kz )  is  the  v  (or  z)  component  of  wave  vector  k, 
and  &0  =  27 t/  A  the  wave  number  in  free  space.  Clearly,  the 
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permittivity  tensor  terms  ( sxx  and  szz)  and  the  permeability 
tensor  term  (juLyy)  are  involved  due  to  the  TM  polarization. 
At  optical  frequencies,  the  permeability  of  natural  materi¬ 
als  can  be  taken  to  unity  [19],  and  therefore  the  effective 
permeability  of  the  ferrofluid  system  is  about  unity.  So  the 
value  of  /uLyy  is  taken  as  unity  in  Eq.  (3)  and  the  following 
analyses.  Negative  refraction  can  be  realized  in  the  ferro¬ 
fluid  system  because  of  the  hyperbolic  equifrequency  sur¬ 
face  and  the  boundary  condition  for  electromagnetic  waves 
at  the  interface  between  a  uniaxial  medium  and  an  iso¬ 
tropic  medium  [26].  One  might  resort  to  Maxwell’s  equa¬ 
tions  to  calculate  the  Poynting  vector  directly.  The  v  and  z 
components  of  the  time-averaged  Poynting  vector  S  are, 
respectively,  given  by  [27] 


d  1.2 


k  H 2 

q,  _  Ax  A1y 

'  V 


ezz  2cos0 


k  H 2 

and  S=^  y 


20JSn 


(4) 


Then  the  angles  of  refraction  for  the  wave  vector  k  and 
Poynting  vector  S  are,  respectively,  given  by 


0k  =  tan  u kx/kz), 


(5) 


6S  =  tan  1  (Sj Sz)  =  tan  \kxexx/(kzezz)\  (6) 

For  the  system  with  Re(g  zz)<0  and  Re(  e»)>  0,  there 
exists  positive  refraction  for  k,  because  of  the  hyperbolic 
equifrequency  contour,  the  conservation  of  the  tangential 
wave  vector,  and  the  causality  theorem  [27].  On  the  basis 
of  Eqs.  (3)  and  (6),  the  effective  refractive  index 


r\  =  sirup/  sin0^ 


(V) 


can  be  computed.  When  Re(ezz)  <  0  and  Re(gxx)  >  0, 
Eq.  (6)  [or  Eq.  (7)]  can  yield  Re(#J  <  0  [or  Re(?7)  <  0], 
namely,  negative  refraction  for  the  Poynting  vector,  as 
shown  in  Fig.  1. 

A  3D  finite  element  method  beyond  the  EM  A. — In  the 
above  analyses,  we  use  the  EMA  [Eq.  (1)]  to  calculate  the 
effective  permittivity  of  the  ferrofluid  system  containing 
chains.  As  a  result,  negative  refraction  has  been  revealed. 
Experimental  results  [11,12]  demonstrated  that  there  exist 
columns  with  equal  spacing  in  confined  ferrofluid  films  in 
the  presence  of  an  in-plane  H.  Accordingly,  we  deduce  that 
there  also  exist  such  columns  with  equal  spacing  in  the 
ferrofluid  system  of  our  interest.  This  should  be  a  reason¬ 
able  deduction  because  similar  to  the  confined  ferrofluid 
films,  our  system  possesses  the  same  dipolar  interaction 
between  (locally  displaced)  chains  that  creates  columns.  In 
order  to  show  the  validity  of  the  above  analyses  using  the 
EMA,  we  conduct  3D  finite  element  simulations  for  a 
ferrofluid  system  containing  columns  with  equal  spacing. 
Certainly  if  we  investigate  the  same  system  with  chains 
rather  than  columns,  in  principle  the  present  3D  simula¬ 
tions  would  work  the  same.  However,  they  are  quite  time- 
consuming.  For  the  comparison  with  Fig.  1  (2D  simula¬ 
tions),  in  Fig.  2  we  also  take  s2  =  1.77  and  the  volume 
fraction  of  the  columns  p  =  0.22.  Each  column  has  a 
length  of  1.0  yam  and  a  permittivity  of  sx  =  —21.4  + 
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FIG.  2  (color  online).  Cross-sectional  view  of  3D  finite  ele¬ 
ment  simulations  of  the  distribution  of  the  absolute  value  of 
electric  field  components  of  an  incident  TM  wave  with  A  = 
758  nm,  in  a  ferrofluid  system  of  thickness  1  yum  containing 
columns  aggregated  by  chains  of  Fe304  nanoparticles  coated  by 
an  Ag  shell.  For  convenience  we  replace  each  column  with  a 
solid  cylinder  by  assuming  the  solid  cylinder  to  possess  the  same 
electromagnetic  responses  as  the  column.  This  assumption  is 
reasonable,  at  least  to  some  extent,  since  the  incident  wavelength 
is  much  larger  than  the  size  of  coated  nanoparticles.  We  set  the 
cross  sections  of  the  cylinders  to  exist  in  a  hexagonal  lattice 
since  equal  spacing  between  columns  was  experimentally  ob¬ 
served  [11].  The  radius  of  the  cylinder  is  30  nm,  and  the  center- 
to-center  separation  between  two  adjacent  cylinders  is  122  nm. 

0.9/.  Also,  a  TM  Gaussian  beam  with  a  width  of  1.6  yarn 
is  incident  at  (p  =  30° .  Figure  2  shows  the  result  simulated 
by  the  3D  finite  element  method  based  on  the  software 
COMSOL  Multiphysics  3.5.  It  clearly  demonstrates  the  ap¬ 
pearance  of  optical  negative  refraction,  which  echoes  with 
Fig.  1.  The  excellent  agreement  is  expected,  since  the 
EMA  is  valid  considering  the  size  and  spacing  of  nano¬ 
particles  are  much  smaller  than  the  wavelength  A.  The 
local  periodic  structures  of  the  chains  or  columns  will 
not  significantly  influence  the  optical  property,  but  can  be 
incorporated  as  a  high-order  approximation  based  on  the 
Ewald-Kornfeld  formulation  [16]. 

Numerical  results. — The  agreement  between  Fig.  1  and 
2  gives  us  confidence  to  further  investigate  the  optical 
properties  of  the  ferrofluid  system  by  using  the  EMA. 
Figure  3  shows  the  real  part  of  rj  as  a  function  of  n,  p, 
and  cp  at  A  =  758  nm.  It  is  worth  reminding  that  higher  H 
leads  to  larger  n.  Evidently  high  H  together  with  suitable  p 
can  lead  to  negative  refraction,  Re(?7)  <  0.  Namely,  one 
can  achieve  the  transition  from  positive  refraction  to  nega¬ 
tive  refraction  for  the  Poynting  vector  by  tuning  H  and/or 
p.  Further  the  negative  refraction  Re(?7)  <  0  is  shown  to 
exist  within  the  full  range  of  (p  (namely,  all-angle  negative 
refraction).  Figure  4  displays  the  imaginary  part  of  77, 
namely  extinction  coefficient  (which  is  proportional  to 
optical  absorption).  We  find  it  varies  with  respect  to  H , 
p ,  and/or  cp.  In  other  words,  one  is  allowed  to  achieve  low 
optical  absorption  by  choosing  H ,  p ,  and  cp  appropriately. 


034501-3 


PRL  104,  034501  (2010) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
22  JANUARY  2010 


FIG.  3  (color  online).  Real  part  of  the  effective  refractive  index 
Rq(tj)  calculated  according  to  Eq.  (7)  as  a  function  of  n ,  p,  and 
</>  (degree)  at  A  =  758  nm.  Three  cross  sections  correspond  to 
n  =  10,  p  =  0.25,  and  <fi  =  70°,  respectively. 

We  also  confirmed  that  (all-angle)  broadband  negative 
refraction  can  be  achieved  in  such  a  ferrofluid  system  by 
investigating  two  other  incident  wavelengths  A  =  661  nm 
and  706  nm  at  n  =  50  and  p  =  0.22  (no  figures  shown 
herein).  For  the  calculations,  we  adopted  the  experimental 
data:  =  -20.09  +  0.45/  [24]  and  ef  =  4.8  +  2.8/  [23] 

at  A  =  661  nm,  and  sa  =  -23.40  +  0.39/  [24]  and  = 
5.0  +  2.4/  [23]  at  A  =  706  nm.  Clearly,  the  wavelength 
bandwidth  can  be  wider  for  the  system. 

In  summary,  for  the  first  time,  we  have  demonstrated 
that  magneto  controllable  all-angle  broadband  negative 
refraction  at  optical  frequencies  can  be  realized  in  aqueous 
ferrofluids,  which  are  made  of  Fe304  nanoparticles  coated 
by  an  Ag  shell.  In  looking  for  experimental  realization,  the 
most  important  point  is  that  the  ratio  between  magnetic 
energy  and  thermal  energy  must  be  high  enough,  which 
causes  the  formation  of  chains  and/or  columns.  Since  the 
permanent  magnetic  moment  of  Fe304  nanoparticles  is 
about  2.4  X  10 4 pc B  [9]  ( i±B :  the  Bohr  magneton),  we  can 
estimate  the  threshold  magnetic  field  14.8  kA/m  above 
which  the  corresponding  magnetic  energy  can  overwhelm 
the  thermal  energy  1/40  eV  at  room  temperature  so  as  to 
obtain  appreciable  //-induced  formation  of  chains/col¬ 
umns.  Practically  one  might  choose  a  higher  H  together 
with  a  lower  temperature.  On  the  other  hand,  the  system  at 
zero  field  ( H  =  0)  may  be  stabilized  against  agglomeration 
by  coating  nanoparticles  with  long  polymer  molecules 


FIG.  4  (color  online).  Same  as  Fig.  3,  but  for  imaginary  part  of 
the  effective  refractive  index  Im(?7). 


(sterically)  or  decorating  them  with  charged  groups  (elec¬ 
trostatically).  The  magnetocontrollable  negative  refraction 
in  ferrofluids  presented  here  promises  a  new  regime  and  a 
number  of  potential  applications  of  metamaterials.  The 
field-induced-assembly  system  of  nanoparticles  may  real¬ 
ize  a  large  specimen  which  is  difficult  to  fabricate  using 
top-down  techniques.  Moreover,  we  can  engineer  the  spa¬ 
tial  dielectric  constant  by  the  external  magnetic  fields  for 
reconfigurable  optical  devices,  such  as  lenses  [27],  invis¬ 
ible  cloaks  [28],  and  waveguides. 
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